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^ Abstract 
>^ 

I We use fits to recent published CPLEAR data on neutral kaon decays to and Trei/ 

to constrain the CPT-violation parameters appearing in a formulation of the neutral 
(— I kaon system as an open quantum-mechanical system. The obtained upper limits of 

the CPT-violation parameters are approaching the range suggested by certain ideas 
concerning quantum gravity. 
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1 Introduction 

The neutral kaon system provides one of the most sensitive laboratories for testing 
discrete symmetries and quantum mechanics at the microscopic level [1]. It is the only 
place where an experimental violation of CP symmetry has yet been seen, and may be used 
to probe the validity of CPT symmetry. This is a property of quantum field theory which 
follows from locality, causality and Lorentz invariance [2]. Violation of CPT would require 
us to revise one or more of these fundamental principles and it is therefore particularly 
important to search for it. Observations of the neutral kaon system have been used to test 
CPT symmetry, in particular within the context of conventional quantum mechanics via 
upper limits on a possible K°-K° mass difference [3, 4, 5]. 

Some approaches to quantum gravity suggest, according to Hawking [6], that quan- 
tum field theory should be modified in such a way that pure quantum-mechanical states 
evolve into mixed states, which would necessarily entail a violation of CPT [7]. This 
possibility may be analyzed using the formulation of open quantum-mechanical systems 
coupled to an unobserved environment [8]. In this framework, the observed system is 
described by a density matrix p that obeys a modified quantum Liouville equation [8] 

p = i[p,H]+^Hp (1) 

where H is the usual quantum-mechanical Hamiltonian. The extra term oc would 
induce a loss of quantum coherence in the observed system, and hence a violation of 
CPT[8, 9, 10]. Since it is conjectured to arise from quantum-gravitational effects, the 

magnitude oi^H may be at most 0(mKVMpi ^ 2x10"^° GeV), where Mpi = 1.2x10^^ GeV 

_ 1^ 

is the gravitational mass scale obtained from Newton's constant: Mpi = G^^ . An equation 
of the form (1) is supported by one interpretation of string theory [11], but could have 
more general applicability. 
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In the case of the neutral kaon system, if the conservation of energy and strangeness 
are assumed, the open-system equation (1) introduces [8] three CPT-violation parame- 
ters a, (3 and 7, which are distinct from the K°-K° mass and decay width differences 
6m = m^o — w^K''5 = Fj^o — Tkc These parameters characterize the openness of an 
apparently isolated kaon to environmental effects which may be associated with quantum 
fluctuations in the space-time background. If present, they would cause the kaon density 
matrix to become mixed, via the suppression of off-diagonal entries in the K°-K° basis, 
thus modifying the rates of Ks and Kl decays as well as interference patterns. In order 
to maintain the hermiticity of the density matrix p and the reality of the entropy, these 
parameters must obey the conditions [8]: 

a, 7 > and a-y > (2) 

Since these parameters would alter the time dependence of the neutral kaon state, they 
are best determined by studying the time evolution of neutral kaon decays, as measured 
by CPLEAR [12, 13]. Measurements [3, 14] at long lifetimes where the loss of quantum 
coherence is increased can be used as additional constraints. The sensitivity reached is in 
the range of interest to the above-mentioned quantum-gravitational ideas. 

2 The CPLEAR measurement of K° and K° decay rate asymmetries 

The CPLEAR experiment [15] is designed to study CP and T violation in the neutral 
kaon system by measuring time dependent decay rate asymmetries of CP and T conjugate 
processes. Initially pure K° and K° states are produced concurrently via the antiproton 
annihilations (pp)j-est ~^ K°K~7r"'" and (pp)j-est ~^ K°K"'"7r~. The strangeness of the 
neutral kaon is tagged by the charge sign of the accompanying kaon. The experimental 
set-up enables us to measure the neutral kaon decay rates in the region of to 20 ts- 
However, the main statistical significance of our measurements is obtained in the region 
below 10 Ts. The correlated analysis of the decays to 7r"'"7r~ and ireu is crucial for the 
determination of the parameters a, /3, 7 as explained below. 

The analysis of the 7r"'"7r~ and ireu final states are described in [12, 13]. The mea- 
sured decay rates into the 7r"'"7r~ final state are shown in figure 1 separately for K° and 
K°, demonstrating visually the known CP violation effect. The decay rate difference is 
measured using the time dependent asymmetry 

^ %o_+,-(r) + axiVKO^.+.-(r) 

In this asymmetry all acceptances common to K° and K° cancel, reducing systematic 
uncertainties. The normalization factor a of the K°, K° rates is fitted simultaneously with 
the CPT violation parameters a, /3, 7 following the procedure described in [12]. The 
measured rates N^o^^+^-{t) and -/Vko^^+^- (r) have been corrected for regeneration [12]. 

For evaluating the yl27r asymmetry (see Eq. 5 below) a precise knowledge of Am 
and the Ks-Kl decay width difference AF is needed. These two parameters are normally 
calculated assuming the time evolution of states according to conventional quantum me- 
chanics. The decay width difference is dominated by the Ks mean life, which is measured 
at early decay times because of statistical reasons, and is therefore not much affected 
by the modifications of the open quantum-mechanical formalism which accumulate with 
time. This is not the case for the mass difference Am and results of experiments using 
conventional quantum mechanical formalism cannot easily be used. However, by using the 
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(3) 



data of the CPLEAR experiment, the value of Am can be obtained in the framework of 
the open quantum-mechanical formalism, by analysing ^Iato [13]: 

^ . ^ ^ [%o^e-.+ I7(^) + iVKO^e+.-.(r)] - [iVKO^e+.-.(^) + ^Vko ^e-.+I7(r )] 

^ [%o_-.+^(r) + iVKO_+.-.(r)] + [iVKO_+.-.(r) + N^o^,-^,^{r)] 

Because of the special form of ^lAm, the systematic errors are reduced and all terms linear 
in the regeneration amplitude cancel [13]. 

3 Description of the Fit 

We use the formalism for the time-dependent K° and K° decay rates to two pions 
and semileptonic final states presented in [16], where the open-system evolution equation 
(1) is solved perturbatively in the small parameters a, /3, 7 and e. Direct CP violation in 
kaon decay amplitudes at the level currently allowed by experiment does not affect our 
analysis and is ignored in the formalism. Since A2Tr{T) contains an important term oc |ep, 
the following full second-order formula for A2Tr{T) is required: 

A2Tr{T) = |2|e| COS ^ + 4/3 sin ^ COS ^ — 8S sin ^ COS ^(|e| sin ^ — 2/3 cos^ ^) 

2a 

cosfAmr — (j) — S(j)) -\ -Xo 

tan (p 



-2^/\e\^ + 0^ cos2(/) ei(i/"^-i/"^)" 
/{l + e('/"^-'/"^)"[7+ |e|' -4;S'cos'(/)-4;S|e|sin(/)]} (5) 



where S, /3, 7 are scaled variables a/AF, /^/AF, 7/AF and the CP impurity parameter e 
is given with the proper phase convention by : 

, , ImMi2 . , , 2Am 

e = e e ^ = , , . : — with <b = atan— — — 

' ' ||AF| +iAm ^ AF 

Furthermore, to the desired order we may write 

2;Scos(^)\ 



8(j) = atan 



Xa = COS 6 (f) sm[AmT — (f)) |AF|rtan^cos(Amr — ^ — 

+ sin (j) cos( Amr — 2(j) — S(j)) 

On the other hand, the following expansion of ylAm(''") to first order, where we assumed 
the validity of the A^ = AQ rule, is sufficient for our purposes: 

2e-|(i/'^5+i/Ti,)T r g _^ 7%(sin Amr - Amr cos Amr)] 

" e-/-^(l + 27) + e-/-s(l - 27) 

We see from equations (5) and (6) that the open quantum mechanical formalism in- 
duces a time dependence which can be distinguished from that of conventional quantum 
mechanics, even with a CPT-violating K°-K° mass difference [5] incorporated. Fitting 
simultaneously the two asymmetries, yl27r (Eq. 5) and ^Iato (Eq. 6) to the CPLEAR data 
(Fig. 2 and Fig. 3), we obtain: 

a = (-0.3 ±3.3) X IQ-^^GeV 
/3 = (2.0 ±2.5) X IQ-'^GeV 
7 = (-0.5± 1.4) X lO-'^GeV 
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|e| = (2.34 ± 0.08) X 10"^ and Am = (526.3 ± 4.2) x lO^/is"^ (7) 

Since the effects of the open quantum-mechanical formalism accumulate with time, the 
best limits on the CPT violation parameters a, /3, 7 are obtained by using in addition 
to the CPLEAR asymmetries, other measurements of Kl decays which refer essentially 

to long decay times [14, 17]. In the formalism used here the parameters \r]^ | and Sl as 

defined in [3] are given for long decay times by: 



2 cos(^ — 28(j)) 



\r,+-r = 7 + \e\' (8) 

cos (p COS"' 0(p 



Sl = 2\e\ ■ COS (j) — 4:f3 COS (j) sin (j) (9) 

Equations (8) and (9) are used as additional constraints, with the values \r]^ | = (2.30 it 

0.035) X 10"^ from [14] and 8l = (3.27 ± 0.12) x 10"^ (average calculated by the Particle 
Data Group [3]). Fitting simultaneously the two CPLEAR asymmetries with the two 
constraints (Eqs. 8 and 9), we obtain: 

a = (-0.5 ±2.8) X 10-^^GeV 
/3 = (2.5 ±2.3) X lO-'^GeV 
7 = (1.1 ±2.5) X lO-^^GeV 



|e| = (2.32 ± 0.06) X 10"^ and Am = (526.3 ± 3.5) x lO^/is"^ (10) 

The main constraint on the parameter 7 is given by the comparison of the measurement 
of the amplitude of the interference pattern (Eq. 5), which defines the value of |e|, with 

the measurement of \r]^ f (Eq. 8), where 7 enters linearly. The result of the fit is shown 

in Fig. 2 and Fig. 3 as a solid line. 

Systematic errors are derived by varying the resolution, the background, and regen- 
eration corrections within their uncertainties [12, 13] and are added in quadrature to the 
statistical errors. The possibility that the observed CP violation is mimicked by the dis- 
cussed CPT violation effects (e = 0) is excluded by the measurements (see also [9, 10, 16]). 

4 Final Results 

The likelihood distribution of the three parameters a, f3 and 7 is approximated by a 
gaussian which takes also into account their correlation, as given by the correlation matrix 
of the fit (see table 1). The 90% confidence level upper limits on each of the parameters 
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Table 1: The correlation coefficients for the parameters a, /3, 7, |e| and Am of the fit 
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are obtained by implementing the positivity constraints (Eq. 2) and integrating over the 
likelihood distribution of the other parameters: 

a < 4.0 X lO-'^GeV < 2.3 x 10"'^ GeV 7 < 3.7 x 10"'' GeV . (11) 

We note that these limits approach the range (9(mK^/Mpi 2 X 10~^° GeV), which is of 
interest for testing ideas based on quantum gravity. 
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2 4 6 8 10 12 14 16 18 20 

Figure 2: Asymmetry yl27r as a function of the decay eigentime. The dots represent the 
CPLEAR data on 7r~'~7r~ decays. The inset displays the data at short decay times. The 
solid line is the result of our fit. The dashed line represents yl27r with positive values of a, 
/3, 7 which are for illustration 10 times larger than the actual limits obtained (Eq. 11). 
The shape of the curve is insensitive to the value of the parameter 7. 
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Figure 3: Asymmetry ^Iato showing the CPLEAR data (dots) on Trez/ decays. The solid 
line is the result of our fit. The dashed line represents ^Iato with positive values of a, /3, 
7 which are for illustration 10 times larger than the actual limits obtained (Eq. 11). The 
shape of the curve depends mainly on the parameter a, apart from Am. 
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